Effects of composition and chemical disorder on the magnetocrystalhne anisotropy of 

Fea,Pti_^. alloys 
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We perform first principles calculations of the magnetocrystalline anisotropy energy (MAE) of the 
Llo-like Fea;Pti_a; samples studied experimentally by Barmak and co-workers in [J. Appl. Phys. 98 
(2005) 033904]. The variation of composition and long-range chemical order in the samples was 
studied in terms of the coherent potential approximation. In accordance with experimental observa- 
tions, we find that, in the presence of long-range chemical disorder, Fe-rich samples exhibit a larger 
MAE than stoichiometric FePt. By considering the site- and species-resolved contributions to the 
MAE, we infer that the MAE is primarily a function of the degree of completeness of the nominal 
Fe layers in the Llo FePt structure. 



Due to its extraordinarily high magnetocrystalline 
anisotropy energy (MAE), LIq FePt is of considerable 
interest to the development of ultrahigh density magnetic 
recording applications, in particular, for heat-assisted 
magnetic recording (HAMR). The LIq phase of FcsoPtso 
is a layered face-centered tetragonal structure, exhibiting 
alternating Fe and Ft layers along the (001) direction. 
FePt also exhibits stable FePta and FegPt phases as well 
as a chemically disordered, cubic phase.'^S Accordingly, 
FePt exhibits phase transitions with respect to compo- 
sition as well as to chemical order and understanding 
the related effects on the magnetic properties is an 
important issue. The large effect of chemical disorder 
on the MAE of FcsoPtso has already been outlined both 
experimentallj0EI and theoretically.'^ 

The degree of long-range chemical order is quanti- 
fied in terms of a chemical order parameter .^El The 
LIq Fe^^Pti-a; alloy is modeled by a repeating sequence 
of two atomic layers, characterized by compositions 
Fcrj^^Pti-rp^ and Fei_rptPtrpt, respectively. The frac- 
tions, rpe and rpt, are related to each other through the 
condition, 1 -I- rpe — Tpt — 2x. Furthermore, set by the 
requirement, rj?e > 1 — rpt (the case of rp^ < 1 — rpt 
can simply be obtained by interchanging the two types 
of layers), the range of rpe is confined to rp,, > x 
{rpt > 1 — x), whereby obviously rpe < min(l,2a;) 
(rpt < min(l,2 — 2a;)). The chemical order parameter 
s is then defined by 



s = 2{rpe -x) ^ 2{rpt ~ I + x) 



(1) 



and ranges from to max(2 — 2a;, 2x). Denoting the 
compositions of the two repeating layers as {A,B), the 
case of complete disorder refers to the compositions 
(Fea;Pti_a;, Fea;Pti_2.) and the maximum order to (Fe, 
Fe2a;-iPt2-22;) for X > 0.5 and to (Fe2xPti-22;,Pt) for 
X < 0.5. Note that only in case of a; = 0.5 can the order 
parameter reach the value s = 1. In the following, we 
refer to the two layers as the nominal Fe layer and the 
nominal Pt layer, respectively. 



TABLE I. Summary of the experimental data obtained for 
four samples of FeajPti-^; in Rcf. 5 

Sample x (%) a (A) c (A) c/a s K (meV/atom) 

1 46.2 3.870 3.721 0.961 0.89 0.453 

2 51.1 3.863 3.710 0.960 0.93 0.709 

3 52.0 3.857 3.706 0.961 0.89 0.775 

4 55.4 3.839 3.704 0.965 0.72 N/A 



Our present study was motivated by the work of 
Barmak and co-workers,'^ who investigated the MAE 
of four Fea:Pti_3^ samples differing in composition and 
degree of chemical order. Table [l] summarizes the 
experimental geometrical and compositional data, as 
well as the measured MAE values for the FePt samples 
studied in Ref. "5", Note that for sample no. 4 the MAE 
could not be determined. One of the main conclusions of 
Ref. [5] is that slightly Fe-rich samples may be preferable 
to FcsoPtso for obtaining a large MAE. In terms of first 
principles calculations we aim to explore the origin of 
this observation, in particular, whether it is it a pure 
effect of composition or whether it is also related to the 
chemical order of the sample. 

To this end, we perform fully relativistic first- 
principles calculations by means of the screened 
Korringa-Kohn-Rostoker (SKKR) method. As the 
method is well documented elsewhere in the literature, 
see e.g. Refs. ITDin^ here we describe only the features 
particularly relevant to this work. We use the local 
spin density approximation (LSDA) of the density 
functional theory (DFT) as parameterized by Vosko et 
s\\^ and treat the potentials within the atomic sphere 
approximation (ASA) . In line with previous work,!-^ the 
self-consistent potentials and fields are calculated from 
scalar-relativistic calculations and the fully relativistic 
Kohn-Sham-Dirac equation is then solved to derive 
the MAE of the system. In all calculations, an overall 
angular-momentum cut-off of ^max = 3 was used. 



The chemical disorder according to the model as 
described above was treated i n term s of the coherent 
potential approximation (CPA)flSIlSl Jt should be men- 
tioned that the MAE of Fea;Coi_a; alloys was recently 
investigated b y usi ng the same model of long-range 
chemical orderP^Ell in the spirit of the magnetic force 
theorem, the MAE is evaluated as the difference in 
the band energy of the system when polarized along 
the easy axis (001) and perpendicular to the easy axis, 
along (100). As in Ref. |7| we estimated the effect of 
temperature-induced spin-fluctuations by scaling down 
the MAE by a factor of 0.6 according to the Langevin 
dynamics simulations of Mryasov et al^ 

In order to verify our method against the experiments, 
first we attempt a direct comparison of our SKKR-CPA 
calculations to the experimental data of Ref. [51 see 
Table |l] As shown in Fig. [l] we performed three sets 
of calculations. The first set only takes into account 
changes in the lattice geometry (i.e., the variation in 
the lattice parameters), while assuming stoichiometric 
composition, x ~ 0.5, and maximum long-range chem- 
ical order, s = 1. Even by taking into account the 
'temperature factor' of 0.6, these calculations yield quite 
high MAE values. Such magnitude differences with 
respect to the experiment are, however, in agreement 
with previous first-principles calculations of the MAE of 
FePt, see e.g. Refs. [20] and [2T1 Furthermore, in this set 
of calculations only a very moderate change (< 3%) of 
the MAE is obtained across the samples. 

In the second set of calculations, we introduce the 
composition x as given in experiment, while keeping the 
degree of chemical order constant at s = 0.89. This 
greatly improves the trend of the MAE, however, the 
relative change of the MAE from sample no. 1 to sample 
no. 2 is still underestimated (< 15%) as compared to 
the experiment (^ 50%). The overall magnitude of 
the MAE is significantly decreased, but it is still by a 
factor of 2-2.5 larger than the measured one. Finally, 
including also the variation of chemical order as given 
in the experiment clearly improves the above mentioned 
relative change between samples no 2 and 3 (~ 30%), 
but, opposite to the experiment, it predicts a slightly 
decreasing trend from sample no. 2 to 3. Note, however, 
that these latter changes are within the range of both 
theoretical and experimental errors. 

Having confirmed that the SKKR-CPA calculations 
satisfactorily reproduce the experimental trends, we next 
consider the general effects of the chemical composition 
X and order parameter s on the MAE of FePt. To this 
end, we used the lattice parameters measured for sample 
no. 3 in Ref. a = 3.857 A and c = 3.706 A, while we 
independently varied the chemical order parameter s as 
well as the composition x. Note that for this theoretical 
study we did not scale down the MAE to mimic temper- 
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FIG. 1. A comparison of the experimental MAE values of 
the Fej;Pti_a: samples studied by Barmak et al. in Ref. 
(open circles) and theoretical MAE values calculated using 
SKKR-CPA as follows, + : using the experimental lattice 
parameters for each sample, but assuming x — 0.5 and s = 1, 
• : using the lattice parameters and the compositions x as 
given in the experiment, but keeping s constant at 0.89, and 
X : using the lattice parameters as well as the values of x and 
s as given in the experiment. Solid lines serve as guide for the 
eyes. 



ature induced effects. The results are shown in Fig. [2] for 
the range of compositions 0.4 < a; < 0.6. (Beyond this 
range, the LIq structure becomes unstable with respect 
to other phases. ^'^) Our results are in good agreement 
with the conclusion of Barmak et aP, inasmuch for any 
given degree of chemical order s the MAE increases 
monotonically with the Fe-content. However, even 
maximally ordered Fej.Pti_j. alloys with a; > 50 cannot 
achieve the MAE of fully ordered FesoPtso (3.31 meV 
per formula unit). 



It should be noted that, at s = 0, the MAE becomes 
negative. This is in contrast to Ref. 22, which reports 
a vanishing MAE for completely disordered FePt under 
the assumption of a cubic unit cell. The 'residual' neg- 
ative MAE we obtain in the case of complete chemical 
disorder is, therefore, due to the tetragonality of the 
lattice (a ^ c). For real samples, where the lattice 
parameters cannot be frozen while varying the chemical 
order and composition, in the case of complete chemical 
disorder the unit cell is expected to become cubic, 
removing thus this 'residual' MAE. 



In order to elucidate the origin of the variation in the 
MAE with the composition and the chemical disorder, 
we consider next the species-resolved contributions to the 
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FIG. 2. The variation of the MAE of Fe^iPti-j; alloys as a 
function of the chemical order parameter s and the composi- 
tion X. Solid lines serve as guide for the eyes. 



MAE. The MAE per unit cell can be decomposed as 



K = rp.D^pt + (1 - rFe)Df,^ 



Pt 



(1 - rpt)DP 



rptD^l 



(2) 

where D^ (/3, 7 = Fe or Pt) denotes the MAE contri- 
bution from an atom of species 7 when it is positioned 
in a nominal /? layer, i.e., within a layer, which in a 
perfectly ordered FesoPtsp alloy would contain only 
atoms of species /3. For the cases of x = 0.4, 0.5 and 0.6, 
in Fig. p^we show D^ as a function of the chemical order 
paramCTer, s. In completely disordered FePt (s = 0), 
the nominal Fe layers and the nominal Pt layers are 
identical. Therefore, at s = 0, the Fe contributions 
in both layers are equal and take a small negative 
value for s — 0, which decreases in magnitude with 
increasing x and practically vanishes at x = 0.6. The 
Pt contributions, on the other hand, are nearly zero for 
all compositions x when s — 0. As the chemical order s 
increases, the Fe contribution in the nominal Fe layers 
rapidly increases up to about 1.8 meV, 2.0 nieV and 2.2 
meV at s = 0.8 for x = 0.4, 0.5 and 0.5, respectively. For 
the fully ordered case, a; = 0.5 and s = I, Dp^ even takes 
the value of about 3.15 meV, close to the total value of 
the MAE (3.31 meV). In contrast, the Fe contribution in 
nominal Pt layers decreases up to s ~ 0.3, then slightly 
increases and, for x > 0.5, reaches a small positive value 
(< 0.5 meV) at maximal chemical order. Remarkably, 
the magnitude of the Pt contributions remain almost 
negligible (< 0.15 meV) over the whole range of chemical 
order. 

As indicated in Fig. |3j the dominant contribution 
to the MAE is rpeD^l, see Eq. ^. It is, therefore, 
intuitive to replot Fig. [2] as a function of the fraction, 
rpf, . This interpretation of the MAE is shown in Fig. |4] 
for different compositions x. Since rpe = x -\- %, it is 
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FIG. 3. Variation in the species-resolved MAE contributions 
against the chemical order parameter a for compositions x = 
0.40 (upper panel), x = 0.50 (middle panel) and x = 0.60 
(lower panel). -I- : Dp"^, contribution of an Fe atom in nominal 
Fe layers, * : Dpf , contribution of a Pt atom in nominal Fe 
layers, x : Dp],, contribution of an Fe atom in nominal Pt 
layers, and 
layers. Solid lines serve as guide for the eyes. 



clear that the horizontal range of the curves in Fig. [2] 
is halved and, more importantly, they are shifted to 
the right by x. As a consequence, for a fixed value 
of rpe the order of the curves with respect to x is 
reversed as compared the order of curves at a given s in 
Fig. [2] . This opposite tendency becomes obvious when 
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FIG. 4. The variation of the MAE of Fe^iPti-j; aUoys as a 
function of the Fe concentration in the nominal Fe layers rpe 
and the overall Fe concentration x. Solid lines serve as guide 
for the eyes. 



considering e.g. the case of rpe = 1 - Fe-rich Y^xPti^x 
will exhibit completely Fe-filled nominal Fe layers at a 
smaller s than Fe5oPt5o, which requires s = 1 in order 
to exhibit completely filled Fe layers. On the other 
hand, increasing disorder (i.e., decreasing s) drastically 



reduces the Fe contribution to the MAE in the nominal 
Fe layer, Dpi, and, consequently, the MAE of the system. 

In conclusion, our calculations strongly support the 
conclusion of Barmak and co-workers in Ref. ,5^ showing 
that, for a given degree of chemical order, the MAE 
increases with the Fe concentration of Fe^^Pti-a;, at 
least within the range 0.4 < x < 0.6. This is due to 
the strongly positive effect on the MAE of the degree of 
Fe- filling of the nominal Fe layers, rpe, which dominates 
the variation in the MAE when varying the composition 
X, while keeping the chemical chemical disorder s con- 
stant. However, Fe^^Pti-a; with x ^ 0.50 cannot attain 
perfect chemical order (s = 1) and perfectly ordered 
FesoPtso yields a larger MAE than the Fe-rich alloys 
with maximum degree of long-range chemical order. 
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